Low-lying excitations and thermodynamics of an antiferromagnetic Heisenberg fractal 

system of a dimension between one and two 



Andreas Voigt *, Johannes Richter and Piotr Tomczak 1 
Institut fur Theoretische Physik 
Otto-von-Guericke-Universitat Magdeburg, Postfach 4120, 39106 Magdeburg, Germany 
Institute of Molecular Physics, Polish Academy of Sciences, 
Smulochowskiego 17, PL-60-179 Poznan, Poland 
(February 1, 2008) 

> , 

£^ • We investigate a frustrated Heisenberg spin-i antiferromagnet on a fractal lattice of dimension 

' d=ln3/ln2 (Sierpiriski gasket). Calculations were performed using (a) exact diagonalization of all 

, eigenstates and eigenvectors for systems up to N=15 and (b) the Decoupled-Cell Quantum- Monte- 

^v^j 1 Carlo method for systems up to N = 366. We present the low-lying spectrum and the specific heat. 

The specific heat shows a second maximum in the low-temperature region. This behavior is similar 
i i . to the behavior of the quantum Heisenberg antiferromagnet on a kagome lattice and suggests a 

' disordered ground state and a spin gap in the considered system. 

^ , Keywords: Heisenberg antiferromagnet, Fractals, specific heat - low temperature 

Introduction. The question for the ground state ordering and the low-temperature thermodynamics of quantum 
antiferromagnets in one or two dimensions receives much interest. In d = 1 there is no long-range order (LRO) in 
the ground state. For d = 2 this question is not generally answered. Various features may influence the magnetic 
ground state ordering: frustration, low coordination number, topological disorder. For several lattices (e.g. square, 



honeycomb, triangular) there is evidence for a long-range ordered ground state, indicating some preference for an 
ordered ground state in d = 2. Only for the kagome lattice the ground state may be disordered. Recently 

S the question has been discussed 'What happens in a quantum antiferromagnet with a dimension between one and 
two?'. In H two of of us (P.T. and J.R.) investigated the ground state of the spin-^ system on the Sierpihski gasket 
| (fractal dimension d = In 3/ In 2) by exact diagonalization and variational Monte-Carlo method. We found rather 
strong arguments that the spin system on the Sierpihski gasket has no long range order and that its ground state may 
O ' correspond to a spin liquid. 

, Additionally a quantum decimation technique was used in jjj by one of us (P.T.) to investigate finite temperature 

properties of that system. It was argued, that the Sierpihski gasket has some similarity with the kagome lattice. 
t-H , In this paper we want to analyze the low-lying spectrum and finite temperature properties of the Sierpihski gasket by 
^ ■ alternative approaches, i.e. quantum Monte Carlo and full diagonalization of small systems. 
The model. We consider the following Hamiltonian 

(N : 

o : H = J E c 1 ) 

t> : 

0^ • where Si is the spin-i operator at the site i and the summation < i,j > runs over all nearest neighbors on a Sierpihski 
fractal lattice. (A system with 366 sites is presented in Fig. 1). The strength of the antiferromagnetic interactions 
is J = 1. The classical ground state of this system is exactly known: The spins belonging to three interpenetrating, 
Q 1 equivalent sublattices form a 120° structure. 

For the quantum system we have calculated all eigenvalues and eigenstates for the systems with N = 6 and TV = 15. 
(The next size is N — 42 and goes beyond all computer facilities available at present). To consider larger systems we 
q have, alternatively, used a decoupled-Cell Quantum Monte-Carlo approach. 

Results. In Fig. 2 we present the low-lying spectrum belonging to different quantum numbers S of the total spin S 2 
for the system with N = 15. According to || a symmetry broken Neel like LRO can be realized in the thermodynamic 
limit (N — » oo) if the lowest levels could be described by an effective Hamiltonian H e ff = Eq + (K/N)S 2 , i.e. the 
^ ■ lowest E a (S) should be linearly dependent on S(S + 1). From Fig. 2 it is obvious that for the Sierpihski gasket there 
is a significant deviation from that linear behavior (cf. straight line between Eo(S — 1/2) and Eq(S = 7/2). Notice, 
that this deviation is even much larger than the corresponding deviation for a linear chain of 16 sites. There is a 
second point of interest concerning ground state ordering: For lattices with antiferromagnetic zero-temperature LRO 
the ground state has minimum total spin (i.e. Sgs — for N even or Sqs = 1/2 for N odd) and the first 
excitation has a total spin S = Sgs + 1. In accordance with the kagome lattice |j]|7) we find here (cf. Fig. 2) (i) a 
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twofold non-trivial degenerated ground state and (ii) 10 additional 5 = 1/2 states between the ground state and the 
lowest S — 3/2 state. Having in mind the above discussion and the variational calculation of the spin correlations in 
0j we argue that for the Sierpihski gasket a symmetry broken Neel like ordering is not favored. 

Another interesting feature of the lowest energies of these system is the existence of two gaps between the lowest 
eigenvalues in the subspace S = 1/2. Numerical calculations show that only these two gaps are necessary in order to 
have two additional peaks in the low-temperature region of the specific heat. 

In Fig. 3 we present the specific heat as a function of the temperature kT/J from the exact diagonalization data for 
N = 6 and N = 15 and from decoupled-cell Quantum-Monte-Carlo calculations for N = 366. This Monte-Carlo 
approach was introduced by Homma et al || and later applied successfully to several quantum spin systems (see 
e.g. ||). For this method the full diagonalization of small clusters is necessary in order to calculate the transition 
probability in the Monte-Carlo process. For the calculations we have used clusters with N = 6, 15. In Fig. 3 we present 
the data with the N = 15 cluster. We find one additional peak in the low temperature region for the system with 
N = 6 and as mentioned above even two such peaks for the system with N = 15. The Quantum-Monte-Carlo data 
for the system with N = 366 also show one additional peak in the low temperature region (The symbols are as high 
as the statistical error of the data). This behavior in the specific heat corresponds to the results obtained by the 
quantum decimation technique and is quite similar to that of a kagome lattice ]To| and supports the suggestion of a 
disordered ground state. 

Conclusion. In general, we argue that the similarity to the kagome lattice may be connected with the lack of LRO 
in both systems. If magnetic correlations are short-ranged the local lattice structure plays the important role. The 
local environment of both lattices is quite similar (four nearest neighbors, triangle as an elementary unit). 
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Fig.l: The Sierpihski gasket with N — 366. 

Fig. 2: The lowest energy eigenvalues of the Sierpihski gasket with N = 15 versus S(S + 1) (S is the quantum number 
of the total spin). 

Fig. 3: The specific heat (per spin) versus temperature. 
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Fig. 1: Voigt et al, Low-lying excitations 




Fig. 2, A.Voigt et a/., Low-lying excitations 



